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In the framework of nonlocal SU(2) chiral quark models with Polyakov loop, we analyze the 
dependence of the deconfinement and chiral restoration critical temperatures on the explicit 
chiral symmetry breaking driven by the current quark mass. Our results are compared with 
those obtained within the standard local Polyakov— Nambu—Jona-Lasinio (PNJL) model 
and with lattice QCD calculations. For a wide range of pion masses, it is found that both 
deconfinement and chiral restoration critical temperatures turn out to be strongly entangled, 
in contrast with the corresponding results within the PNJL model. In addition, it is seen 
that the growth of the critical temperatures with the pion mass above the physical point is 
basically linear, with a slope parameter which is close to the existing lattice QCD estimates. 
On the other hand, within the present mean field calculation we find an early onset of the 
first order transition expected in the large quark mass limit. 

PACS numbers: 12.39.Ki, 11.30.Rd, 12.38.Mh 



I. INTRODUCTION 



It is widely believed that as the temperature and/or density increase, strongly interacting matter 
undergoes some kind of transition from a hadronic phase, in which chiral symmetry is broken and 
quarks are confined, to a partonic phase in which chiral symmetry is restored and/or quarks are 
deconfined. The detailed understanding of this phenomenon is relevant not only in particle physics 
but also e.g. in the the study of the early universe, the interior of neutron stars, etc., therefore it has 
become an issue of great interest in recent years, both theoretically and experimentally From the 
theoretical point of view, one way to address this problem is through lattice QCD calculations 
J]. However, even if significant improvements have been done in this field in the last years, this 
ab initio approach is not yet able to provide a full understanding of the QCD phase diagram. One 
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serious difficulty in this sense is given by the so-called sign problem, which prevents straightforward 
simulations at finite baryon density. In this situation it is worth to develop alternative approaches, 
such as the study of effective models that show consistency with lattice QCD results and can 
be extrapolated into regions not accessible by lattice techniques. Here we will concentrate on 
one particular class of effective theories, namely the so-called nonlocal Polyakov— Nambu— Jona- 
Lasinio (nlPNJL) models J 8], in which quarks move in a background color field and interact 
through covariant nonlocal chirally symmetric four point couplings. Related Polyakov-Dyson- 



■ |£|. These approaches, which can be 
ICHla]. offer a common framework to 



Schwinger equation models have also been recently analyzec 
considered as an improvement over the (local) PNJL model 
study both the chiral restoration and deconfinement transitions. In fact, the nonlocal character 
of the interactions arises naturally in the context of several successful approaches to low-energy 
, and leads to a momentum dependence in the quark propagator that can 



quark dynamics 
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be made consistent 



19( 1 with lattice results 



20 



211 ]. Moreover, it has been found that, under 



certain conditions, it is possible to derive the main features of nlPNJL models starting directly 

'3-2 



from QCD 



that nonlocal models 



22]. From the phenomenological side, it has been shown 
provide a satisfactory description of hadron properties at zero temperature and density. 

As mentioned above, it is important to consider situations in which the results obtained within 
effective models can be compared with available lattice QCD calculations. For example, it is 
clear that vacuum properties such as the pion mass and decay constant, as well as other features 
related to the chiral/deconfinement transitions (like e.g. the nature of the transitions, or the critical 
temperatures) will depend on basic parameters of QCD, such as the number of quark flavors and the 
values of current quark masses m^. In particular, for the simplified case of two degenerate flavors 
with ruu = rrid = m, the dependence of several relevant quantities on m has been studied with 
some detail in lattice QCD. Thus, the corresponding analysis within nlPNJL models can provide 
an interesting test of the reliability of this effective approach. Actually, it has already been shown 



that several chiral effective models 



27 



291 ] are not able to reproduce the behavior of the critical 



temperatures observed in lattice QCD when one varies the parameters that explicitly break chiral 
symmetry (i.e. the current quark masses, or the pion mass in the case of meson models) at vanishing 
chemical potential. This fact has been taken as an indication that the transition may be not just 
dominated by pure chiral dynamics [s^ . It is worth to notice that in the framework of the standard 
(local) NJL model the enhancement of the critical temperature with m is too strong in comparison 
with lattice QCD estimates. Although the inclusion of confinement effects through the coupling 
to the Polyakov loop weakens this enhancement, one finds a too large splitting between the chiral 
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restoration and deconfinement transition temperatures [31[. The presence of confinement effects 
togetlier witli a strong entanglement between tlie chiral restoration and deconfinement transitions 
is indeed one of the features of nlPN JL models [32] • 



In view of the above mentioned points, the aim of the present work is to study the effect of 
explicit chiral symmetry breaking on the deconfinement and chiral restoration critical temperatures 
within nlPNJL models. This article is organized as follows. In Sec. II we provide a description of 
the model, proposing two alternative parameterizations. In Sec. Ill we analyze the m-dependence 
of some pion properties and compare the results with existing lattice calculations. In Sec. IV 
we analyze the current quark mass dependence of the critical temperatures at vanishing chemical 
potential, comparing our results with those obtained in alternative models and lattice QCD. Finally 
in Sec. V we summarize our main results and conclusions. 

II. FORMALISM 

We consider a nonlocal SU(2) chiral quark model that includes quark couplings to the color 



gauge fields. The corresponding Euclidean effective action is given by 



33| 



Se = I d^x I V^(x) {-ii^^D^ + m) ^{x) - % ja{x)ja{x) - jp{x)jp{x) + U \ , (1) 
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where is the = 2 fermion doublet ^ = {u,d)'^, and rh = diag{mu,md) is the current quark 
mass matrix. In what follows we consider isospin symmetry, = = m. The fermion kinetic 
term in Eq. ([1]) includes a covariant derivative = — iA^^ where are color gauge fields, and 
the operator 7^9^ in Euclidean space is defined as 7 • V + ^^d/dr, with 74 = ^79. The nonlocal 
currents jaix), jp{x) are given by 



jaix) = J d^Z G{z) i^{x + ^ Ta^(x-^ 



,4 ^, s T / z\ i $ , / z 



jp{x) = J d'^z F{z) ^ (^x - -j , (2) 

^ y 

where, Ta = {t^i^^r) and u{x') d v{x) = u{x')dxv{x) — dx'u{x')v{x). The functions Q{z) and 
J-{z) in Eq. ([5]) are nonlocal covariant form factors characterizing the corresponding interactions. 
Notice that the four currents ja{x) require a common form factor Q{z) in order to guarantee 
chiral invariance, while the coupling jp{x)jp{x) is self-invariant under chiral transformations. The 
scalar-isoscalar component of the ja{x) current will generate a momentum dependent quark mass 
in the quark propagator, while the "momentum" current jp{x) will be responsible for a momentum 
dependent quark wave function renormalization. Now we perform a bosonization of the theory, 
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introducing bosonic fields a'i^2ix) and TTa{x), and integrating out the quark fields. Details of this 
procedure as well as of the determination of vacuum and meson properties at vanishing temperature 
in this framework can be found e.g. in Ref. [3]. 

Since we are interested in the deconfinement and chiral restoration critical temperatures, we 
extend the bosonized effective action to finite temperature T. This can be done by using the 
standard Matsubara formalism. Concerning the gauge fields A^, we assume that quarks move on 
a constant background field (p = A4 = iAq = ig 6^q Ga A"/2, where Ga are SU(3) color gauge fields. 
Then the traced Polyakov loop, which in the infinite quark mass limit can be taken as an order 
parameter of confinement, is given by $ = |Tr exp(i(/)/T). We work in the so-called Polyakov 
gauge, in which the matrix (p is given a diagonal representation (p = 03A3 + (ps^s- This leaves 
only two independent variables, (ps and 08 . In the case of vanishing chemical potential, owing to 
the charge conjugation properties of the QCD lagrangian, the mean field traced Polyakov loop is 
expected to be a real quantity. Since (p^ and have to be real valued, this condition implies 
4>s = 0. The mean field traced Polyakov loop reads then $ = $* = [1 + 2 cos {(p-^/T)] /3. Thus 
in the mean field approximation, which will be used throughout this work, the thermodynamical 
potential O'^^^ at finite temperature and zero chemical potential is given by 



MFA 



E 

c=T,g,b n=— 00 



log 



(27r)3 

where M(p) and Z{p) are given by 

M{p) = Z{p) [m, + cji g{p)] , 



+ 



^1 + 

2Gs 



Z{p) = [l-a2 f{p)]-' 



+ Z^($,<I>*,r) , (3) 



(4) 



Here (Ti^2 are the mean field values of the scalar fields (note that iTa = 0), while and f{p), g{p) are 
Fourier transforms of J-{z) and Q{z), respectively. We have also defined 



Pn,p 



(2n+ l)7rr + 



+ p' 



(5) 



where the quantities (pc are given by the relation cp = dia.g{(pr, (pg, (ph) = diag((/)3, —(p^,Q). 

To proceed we need to specify the explicit form of the Polyakov loop effective potential 
<!>*, T). We consider two alternative functional forms commonly used in the literature. The 
first one, based on a Ginzburg-Landau ansatz, reads 3] 



where 



|(|)|2 _|_ |$*|2~] 



Q V V 7 y 16 ^' ' 



1$ 



*\2\ 



T 



+ a2 



+ 0.3 



(6) 



(7) 
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The potential parameters can be fitted to pure gauge lattice QCD data so as to properly reproduce 



the corresponding equation of state and Polyakov loop behavior. This yields 13 1 



ao = 6.75 , ai = —1.95 , a2 = 2.625 
a3 = -7.44 , b-s = 0.75 , 64 = 7.5 . 



(8) 



A second usual form is based on the logarithmic expression of the 



iaar measure associated with 



the SU(3) color group integration. The potential reads in this case 1^ 

Uiog{<^, T) = I - ^ a{T) m* + b{T) log [1 - 6 + 4 $3 + 4 ($*)3 _ 3 ($$*)2] | j-^ ^ (9) 
where the coefficients are parameterized as 

a{T) = ao + aJ^)+aJ^y , b{T) = bJ^Y . (10) 



Once again the values of the constants can be fitted to pure gauge lattice QCD results. This leads 
to Q 



ao = 3.51 



ai 



-2.47 



02 = 15.2 



-1.75 



(11) 



The dimensionful parameter Tq in Eqs. ([7]) and (jlOp corresponds in principle to the deconfinement 
transition temperature in the pure Yang-Mills theory, Tq = 270 MeV. However, it has been ar- 
gued that in the presence of light dynamical quarks this temperature scale should be adequately 
reduced [341 ]. 



Finally, one has to take into account that fi'^^^ turns out to be divergent, thus it has to be 
regularized. Here we use the prescription described e.g. in Ref. 



35l |. namely 



n 



UFA 



n 



MFA 



free 



(12) 



where Q^'^'^'^ is obtained from Eq. by setting cti = (T2 = 0, and n^^^g is the regularized expression 
for the quark thermodynamical potential in the absence of the four point fermion interaction, 



free 
reg 



-4r 



(2^ 



1 + exp 



T 



(13) 



c=r,g,b s=±l 

with Cp = 2 _|_ jj^2 ^ "Yhe last term in Eq. ()12p is just a constant fixed by the condition that 



n^l^ vanishes at T = 0. 



Given the full form of the thermodynamical potential, the mean field values ai^2 and (p^ can be 
obtained as solutions of the coupled set of "gap equations" 



oqMFA 



{dai,da2,d(j)3) 







(14) 
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Once these mean field values are obtained, the behavior of other relevant quantities as functions 
of the temperature and chemical potential can be determined. We concentrate in particular in 
the chiral quark condensate {qq) = d^^^^/dm and the traced Polyakov loop which will be 
taken as order parameters of the chiral restoration and deconfinement transitions, respectively. 
The associated susceptibilities will be defined as Xch = d {qq)/dm and xpl = d^/dT. 

In order to fully specify the model under consideration we proceed to fix the model parameters 
as well as the nonlocal form factors g{q) and f{q) at the physical point m.^ = m^^^ = 139 MeV. 
We consider two different functional dependences for the form factors. The first one corresponds 
to the often used exponential functions 

g{q) = exp {-q^/Al) , f{q) = exp (-^VA?) , (15) 

which guarantee a fast ultraviolet convergence of the loop integrals. Note that the range (in 
momentum space) of the nonlocality in each channel is determined by the parameters Aq and 
Ai, respectively. Fixing the current quark mass and chiral quark condensate at T = = 
to the phenomenologically adequate values m = 5.7 MeV and (qq)^^^ = 240 MeV, the rest of the 
parameters can be determined so as to reproduce the physical values of /jr and ttt-tt, and by requiring 



Z{0) = 0.7, which is within the range of values suggested by recent lattice calculations [20|, l2ll ]. 
In what follows this choice of model parameters and form factors will be referred to as SI. The 
second type of form factor functional forms considered here is given by 

where 

fm{q) = [l+{q'/^lfY' ' = [1+ («VA?)]"'^'- (17) 

As shown in Ref. Q, taking m = 2.37 MeV, = 309 MeV, a, = -0.3, Aq = 850 MeV and 
Ai = 1400 MeV one can very well reproduce the momentum dependence of mass and wave function 
renormalization obtained in lattice calculations, as well as the physical values of rriT^ and 7,^. In 
what follows this choice of model parameters and form factors will be referred to as S2. Details on 
the model parameters and the predictions for several meson properties in vacuum can be found in 
Ref. 



f. Q. 
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III. ZERO TEMPERATURE PSEUDOSCALAR MASS AND DECAY CONSTANT 

AWAY FROM THE PHYSICAL POINT 

As stated, we want to study the dependence of nlPNJL model predictions on the amount of 
exphcit chiral symmetry breaking. This can be addressed by varying the current quark mass m, 
while keeping the rest of the model parameters fixed at their values at the physical point. As a first 
step we analyze in this section the corresponding behavior of the pion mass and decay constant at 
vanishing temperature, in comparison with that obtained in the (local) NJL model and in lattice 
QCD. Our results are shown in Fig. 1. As it is usual in lattice QCD literature, we choose to take 
ttItt instead of m as the independent variable in the plots. The main reason for this is that 
is an observable, i.e. a scale independent quantity, whereas m is scale dependent, hence its value 
is subject to possible ambiguities related to the choice of the renormalization point. Dashed and 
solid lines correspond to parameter sets SI and S2, respectively, while dotted lines correspond to 
the curves obtained within the NJL model. Fat dots stand for lattice QCD results from Ref. [36|. 



The upper panel shows the behavior of the ratio m^/m as a function of niT^. In order to account 
for the above mentioned renormalization point ambiguities, the corresponding quark masses have 



been normalized so as to yield the lattice value w-J^^ ~ 4.452 MeV at the physical point (36|]. From 
the figure one observes that both NJL and nlPNJL models reproduce qualitatively the results 
from lattice QCD, showing a particularly good agreement in the case of the nlPNJL model for 
parameter set S2. However, the situation is different in the case of f-,^ (lower panel in Fig. 1): 
while the predictions from nonlocal models follow a steady increase with m-,^, in agreement with 
lattice results, the local NJL model badly fails to reproduce this behavior. Moreover, it can be 
seen that the discrepancy cannot be cured even if one allows the coupling Gs to depend on the 



3l[. Thus, these results 



current quark mass (we have taken Gs as a constant in nlPNJL models) 
can be considered as a further indication in favor of the inclusion of nonlocal interactions as a step 
towards a more realistic description of low momenta QCD dynamics. 

IV. DEPENDENCE OF CRITICAL TEMPERATURES ON EXPLICIT CHIRAL 

SYMMETRY BREAKING 

In this section we analyze within our nonlocal models the mass dependence of the critical tem- 
peratures for the deconfinement and chiral restoration transitions at vanishing chemical potential. 
We start by considering the temperature dependence of the chiral and deconfinement order pa- 
rameters, as well as the corresponding susceptibilities, for some representative values of the pion 
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mass. The corresponding results for the lattice motivated parameterization S2 are shown in Fig. [2l 
including both the case of the polynomic (left panels) and logarithmic (right panels) Polyakov 
potentials. Qualitatively similar results are found for the exponential parameterization SI. Let us 
first discuss the results for the polynomic potential. From the figure it is seen that both transitions 
proceed as smooth crossovers, as expected from lattice QCD results. Moreover, we observe that 
as ttItt increases, the position of the peaks of the susceptibilities Xch and xpl (left lower panel) 
move simultaneously towards higher values of T, the difference between the corresponding critical 
temperatures being in all cases at the level of a few MeV. It is also seen that as increases 
the chiral restoration transition tends to be less pronounced, while the confinement one becomes 
steeper. In the case of the logarithmic potential, we also observe that the transition tempera- 
tures increase with m-,^, as expected. However, for a given value of m-,^ both the chiral restoration 
and deconfinement transitions are steeper than in the case of the polynomic potential, and the 
correlation between them is stronger (e.g. the difference between the transition temperatures for 
= m^^^ is now about 0.02 MeV). In fact, it turns out that already for m-j^ = 500 MeV one 

finds a first order phase transition. From lattice QCD results the onset of a first order phase 

□ 

transition is indeed expected above a certain critical amount of explicit symmetry breaking |37l |. 



However, present estimations [38H4[lt| indicate that the corresponding critical pseudoscalar mass 
should be much larger than the physical pion mass; therefore, the early change in the character 
of the transition appears as an unrealistic feature of the logarithmic Polyakov potential. In the 
case of the polynomic potential the onset of this first order phase transition occurs for a pion mass 
larger than 700 MeV, i.e. in an energy region where the applicability of the effective quark models 
is limited. This situation is qualitatively similar to that observed in the local PNJL model. 

In Fig. [3] we show the results for the mass dependence of the critical transition temperatures 
within our nonlocal models. For comparison we also quote typical curves obtained in the framework 
of the local PNJL model (here we have considered the parameterization in Ref. 3]). Upper 
and lower panels correspond to polynomic and logarithmic Polyakov potentials, respectively, with 
To = 270 MeV. Before discussing in detail the results obtained for the nlPNJL models, let us 
comment those corresponding to the PNJL model: from Fig. [3] we observe that already at the 
physical value = m^^^ the model predicts a noticeable splitting between the chiral restoration 
temperature T^h (dashed line) and the deconfinement temperature TpL (dotted line). In addition, it 
is seen that the growth of Tch with is stronger than that of TpL , which implies that the splitting 
between both critical temperatures becomes larger if m-j^ is increased. This is not supported by 
existing lattice results 
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42( 1 ■ which indicate that both transitions take place at approximately 
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the same temperature, up to values of 171.,^ even larger than those considered here. Comparing 
both panels it is seen that the splitting is more pronounced for the PNJL model that includes a 
logarithmic Polyakov potential. 

We turn now to the curves obtained within nonlocal models. First of all, from the figure it 
is seen that both parameterizations SI and S2 lead to qualitatively similar results. Contrary to 
the situation in the PNJL model, in nlPNJL models both the chiral restoration and deconfinement 
transitions occur at basically the same temperature for all considered values of rriT^. Moreover, com- 
paring the results for the two alternative Polyakov loop potentials we see that the main qualitative 
difference between them is the already mentioned fact that in the case of the logarithmic potential 
there is a critical pion mass of about 400 MeV where the character of the transition changes from 
crossover to first order (dashed-dotted and solid lines in the lower panel of Fig. [3l respectively). 
By analyzing in more detail the pion mass dependence of the critical temperatures, it is seen that 
for m^r above the physical mass the nlPNJL model results can be accurately adjusted through a 
linear function 



This is in agreement with the findings of the lattice calculations of Refs. 



(18) 



41 



42( 1 . Our results 



for the slope parameter A for both parameterizations and Polyakov loop potentials are in the 



range of 0.06 — 0.07. For comparison, most lattice calculations find A < 0.05 
according to some recent analyses 
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43|-l45|], while 



461 ] the value could be somewhat above this bound. Thus 
the slope parameter predicted by the nonlocal PNJL models appears to be compatible with lattice 
estimates. This can be contrasted with the results obtained within pure chiral models, where one 
finds a strong increase of the chiral restoration temperature with m,, 
the chiral quark model of Ref. 



27 



291 ]. For example, within 



291] one gets a value A = 0.243. 
It is also worth to discuss the effect of considering a value of Tq that depends on the current quark 
masses, as suggested in Ref. 



341 ]. The change in Tq leads to an overall decrease of the transition 
temperatures, which keep the rising linear dependence on m,r but with a slope parameter that 
gets reduced by about 15 — 20%. The main noticeable difference is that in all cases the transition 
becomes steeper, which leads to an earlier onset of the first order transition. For example, for the 
parameter set S2 we find that the transition becomes of first order already at ~ 500 MeV in the 
case of the polynomic Polyakov potential, and about one half of this value for the logarithmic one. 
These critical masses appear to be too small in comparison with present lattice QCD estimations. 



However, in this respect it is important to recall the importance of considering corrections that go 
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beyond the mean field approximation used here. In fact, although the role of these corrections is 
expected to be less important as the quark mass increases Q], in the range of masses considered here 
they can be significant enough to soften the transitions and lower the critical temperatures Q]. In 
this sense, although a fully nonperturbative scheme to account for meson fluctuations in nonlocal 



models is still lacking, some important steps have been taken 



Ref. 



As it is pointed out in 



48l |. these fluctuations could also help to avoid thermodynamical instabilities that could arise 



in nonlocal models. 



V. SUMMARY AND CONCLUSIONS 



In this work we have analyzed the dependence of the deconfinement and chiral restoration critical 
temperatures on the explicit chiral symmetry breaking driven by the current quark mass. We work 
in the framework of SU(2) nonlocal chiral quark models with Polyakov loop (nlPNJL models), 
considering two different functional forms of the Polyakov loop effective potential commonly used 
in the literature, namely a polynomic function and a logarithmic function. As a first step we have 
considered the mass dependence of the pion mass and decay constant at vanishing temperature, in 
comparison with that obtained in the local NJL model and in lattice QCD. We have found that, 
while lattice results for the ratio m^/m are in agreement with both local and nonlocal models, 
those for show a significant increase with ttt-tt that can be reproduced only by the predictions 
of nonlocal models. Concerning the deconfinement and chiral restoration critical temperatures, we 
have found that, contrary to the case of the local PNJL model, in nlPNJL models both critical 
temperatures turn out to be strongly entangled for the considered range of pion masses. In addition, 
it is seen that the growth of critical temperatures with the pion mass above the physical point is 
basically linear, with a slope parameter which is close to existing lattice QCD estimates. On the 
other hand, particularly in the case of the logarithmic Polyakov loop potential, the present mean 
field calculation leads to a too early onset of the first order transition known to exist in the large 
quark mass limit. We expect that the development of a fully nonperturbative scheme to account 
for meson fluctuations in nonlocal models might help to provide a solution to this problem. 
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FIG. 1: Pion properties at T = as functions of the pion mass in local and nonlocal chiral quark models. 
Upper an lower panels correspond to the ratio rn^/nic and the pion decay constant /^r, respectively. Lattice 
results are taken from Ref. 
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FIG. 2: Order parameters (upper panels) and the eorresponding susceptibilities (lower panels) as functions 
of the temeprature for some representative values of the pion mass. Left (right) panels correspond to the 
polynomic (logarithmic) Polyakov potential. 
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FIG. 3: Critical temperatures as functions of the pion mass for PNJL and nlPNJL models, considering 
polynomic (upper panel) and logarithmic (lower panel) Polyakov loop potentials. Dashed and dotted lines 
correspond to chiral restoration and deconfinement transition temperatures, respectively. For the nlPNJL 
models with a logarithmic potential (lower panel), both transitions occur at the same temperature, and they 
can be of first order (solid lines) or proceed as a smooth crossover (dashed-dotted lines). 



